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Abstract 

The dynamics of a packages diffusion process within a selforganized 
network is analytically studied by means of an extended /-spin fa- 
cilitated kinetic Ising model (Fredrickson-Andersen model) using a 
Fock-space representation for the master equation. To map the three 
component system (active, passive and packages cells) onto a lattice 
we apply two types of second quantized operators. The active cells 
correspond to mobile states whereas the passive cells correspond to 
immobile states of the Fredrickson-Andersen model. An inherent co- 
operativity is included assuming that the local dynamics and subse- 
quently the local mobilities are restricted by the occupation of neigh- 
boring cells. Depending on a temperature-like parameter h^ 1 (inter- 
connectivity) the diffusive process of the packages (information) can 
be almost stopped, thus we get a well separation of the time regimes 
and a quasi-localization for the intermediate range at low tempera- 
tures. 

1 Introduction 

During the last years there is a huge effort to understand the kinetics of non- 
equilibrium phenomena. A wide range of discrete and continuous models 
for such processes is analytically and numerically studied. The problems of 
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interest in this context concern the crystal growth, transport (traffic) models, 
diffusion processes and supercooled liquids |IJ-§|. 

Here, we will apply the kinetics of the Fredrickson- Andersen model (FAM) 
recently discussed in the framework of the glass transition and related phe- 
nomena PHl2|- But we show that this model may be used on other fields 
like stock trading, citation networks, company relations or internet commu- 



nications as well |13Hl6l. In general, we study the diffusion of information 
within a network system of active links and passive/active cells (or nodes). 
The switch between a passive and active cell is controlled by a temperature- 
like parameter /i _1 which can be interpreted as the interconnectivity. At the 
maximum interconnectivity (at infinite temperature) the system possesses 
equal parts of active and passive cells (nodes) whereas at the minimum in- 
terconnectivity (at zero temperature) there are only passive cells (nodes). 
Furthermore, the alteration is also controlled by the nearest neighborhood. 
If enough adjacent active cells (nodes) exist (more than a fixed number /) 
the active cell (node) can become passive and vice versa. Thus, only a suffi- 
ciently active environment may determine and alter the state of a cell (node) 
like in a citation community where only accepted (active) people may de- 
cide about the worth of an opinion of a member in a related field. To this 
network formed by passive/active cells (nodes) and active links, consisting 
of two adjacent active cells (nodes), we add further particles which may be 
assigned to (information) packages. These particles can only diffuse along 
active links but are confined by passive cells (nodes). Therefore, we have a 
diffusion in a self-organized network where information may stick (and there- 
fore be localized) at passive cells (nodes) or runs through a network of active 
links. This is the same situation like in the internet where passive routers 
cannot transfer any data but data on active router should be directed to the 
next active router (Of course, the difference is that this is a directed motion 
whereas we consider a non-directed motion here.). 

Now we demonstrate how one may relate this network to the FAM origi- 
nally dedicated to supercooled liquids. For such systems one assumes that the 
dynamics can be roughly regarded as independent of the microscopic details, 
thus the analysis is based on a mesoscopic formulation. Hence, one reduces 
the degrees of freedom to a smaller set of relevant observables following the 
original idea of Fredrickson and Andersen J?]-[|. The supercooled liquid as 
the many-body system is treated on a virtual lattice with sufficiently small 
cells of size I where the lattice should not have any significant influence on 
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the underlying dynamics. However, in case of our network model we will 
immediately start on the mesoscopic scale (e.g. on the elementary scale of 
routers where the macroscopic scale is the total internet). 

The cell structure enables us to attach to each cell j a two-state observ- 
able CTj which characterizes its actual state. The realization is given by the 
local activity rrij with Oj = —1 (spin down) if rrij < fh (passive state) and 
<jj = +1 (spin up) if rrij > fh (active state) where m is a fixed threshold 
of the system. The threshold is chosen in this way that all cells are passive 
at zero temperature. Active (passive) states may be also denoted as mobile 
(immobile) states. The number of cells in the active and the passive state, 
respectively, have not to be conserved. (Notice that in case of supercooled 
liquids the usual separation is done by two different density states although 
this is more connected to a Kawasaki exchange dynamics rather than to a 
Glauber flip dynamics |J. This mapping implies different densities of such 
cells; Oj = — 1 corresponds to the more solid-like state and Oj = 1 corresponds 
to the more liquid-like state.) Thus, the real kinetics is more Glauber-like 
(non-conserved) dynamics thus we suppose that the basic dynamics is a flip 
process o~j = +1 <->• Oj = — 1. It is controlled by self-induced topological 
restrictions introduced below. This type of dynamics leads to a relaxation 
behavior resembled to that of a cooperative system. In particular, an elemen- 
tary switch at a given cell is allowed only if the number of nearest neighbor 
mobile/active cells (oj = +1) is equal to or larger than a certain restrictive 
number / with < / < z (z is the coordination number of the lattice). 
Hence, elementary flip processes combined with the restrictions may lead to 
the cooperative arrangements within the underlying network system. 

As demonstrated such a mesoscopic approach is able to model supercooled 
liquids as well. The Fredrickson- Andersen model (or /-spin facilitated kinetic 
Ising model) [0-0] was analytically P,|6|,[i7|| and numerically p^ HT2] , [T8|] studied 
in broad variety. The FAM can be classified as an Ising-like model of which 
kinetics is confined by restrictions of the ordering of nearest neighbors to 
a given lattice cell. These self-adapting environments especially influence 
the long-time behavior of the spin-relaxation |B|,[TT| particularly relevant for 
supercooled liquids. 

Additionally, here we add particles to the system of which diffusive dy- 
namics are coupled with the existence of the active states. Therefore, we 
require that active cells are necessary for the motion of packages because 
passive states block any packages diffusion. Or in other words data are fixed 
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and static as long as they are assigned to passive cells, but information may 
diffuse through a network of active cells. In the present work we will incorpo- 
rate this features into the Fock-space representation of the master equation in 
order to compute the packages concentration in a continuous mean-field- like 
approximation. 

2 Algebraic Representation and Master Equa- 
tion 

Firstly, we give a short survey about the Fock-space method (known as quan- 
tum Hamiltonian as well - more details can be found in e.g. ||) A certain 
state in our lattice system can be characterized by a set of discrete numbers 
n = {rii} (or resp. v = {vi}) where n^Vi G {0, 1} denote the local state of 
a lattice cell i. Furthermore, the following convention is used: rii — 1 (0) is 
related to the active (passive) state, also denoted by A (B). The state v j = 1 
(0) corresponds to a cell % occupied (non-occupied) by a package. We start 
from the one-step master equation formulated in the following form 

d t P(n,v,t) = L'P{n,v,t) (1) 

where P(n,v,t) is the probability for a certain configuration {n, v}. The 
linear operator L' specified by the dynamics of the system describes the time 
evolution. Then, following |19|,P0|, the probabilities P(n,v,t) can be related 
to the Fock-space state vector | F(t)) as a weight using the decomposition 
into the base vectors | n)<g> | v) of a orthonormal vector space, 

I F(t)) = Y,P(n,v,t) |n>® | v). (2) 

n,v 

This equation enables the quantum formulation for the master equation re- 
sulting in an equivalent manner 

d t | F(t)) = L | F(t)} (3) 

where L' is assigned to the operator L in the Fock-space representation. The 
procedure had been originally derived for Bose-like system and was 



later applied also to Fermi-like systems ||,|2T],^2|]. Recently, we have proposed 
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a further extension applying Para-Fermi statistics for different restricted oc- 
cupation numbers |||!|]. The average of a physical quantity G(n,v) is given 
by the trace over G: 

(G(t)) = J2 p (n, v, t)G{n, v) = (r\G\ F(t)) (4) 

n,v 

where (r | is the reference state related to the base vectors by 

<t ? |=X>I®<*I=®( \ )■ (5) 

ri.v \ / 

The reference state is completely determined by the base { | n),\ v)} of the 
Fock space and not assigned to the special model or to the evolution operator 
L. The conservation of the total probability is manifested by (r | L = 0. 
Therefore, the equation of motion is given by 

F(t)). (6) 

Notice that the dynamical equation depends on both the algebraic properties 
of the underlying operators and the mathematical structure of L. Thus, we 
introduce the second quantized lowering Oi(vi) and raising a\(v\) operators 
forming the evolution operator L to create the base states | n) (| v)) from a 
vacuum state | 0). Both types of (independent from each other and hence 
commuting) operators fulfil the relationship (for Paulions) 

a^a] + ajdj = 8^ + 2a]aj (1 — 8^) . (7) 

As mentioned above, the inherent properties of the FAM is the restriction of 
the flip dynamics at cell i, Oi <-> — a,- n 

2 EK' I (! + Oj) I n o) = E(% I Ai I n j) > f (8) 

where j(z) means a sum over all adjacent cells of i and / is the restriction 
number. The number operators Aj and Vj denote OjCij and VjVj as usual. 
Concerning the motion of information packages we postulate a diffusive mo- 
tion of the particles coupled to the existence of active cells at both related 



d t (G(t)) = (f\GL\ F(t)) = (f | G,L 
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sites. This exchange process due to Kawasaki enhances the mobility of pack- 
ages in active neighborhoods whereas it slows down inside a passive cluster. 
Summarizing we obtain the evolution operator taking the form 

L = L F + L E (9) 
Lf = +^2\ba(1 - ai)a\ ^ A mi ■ ■ ■ A mf 

i <m\---mf,i> 
i <rri\---mf ,i> 

L E = + D Q v\.v s — (l — V r ) V s A r A s + symmetric term 



<rs> 



with the quantities \ab = Aexp[/i], Aba = Aexp[— h] and D are the ki- 
netic coefficients of the diffusion process. They are appropriate thermody- 
namic weighted to fulfil the detailed balance condition. The parameter h 
corresponds to the inverse temperature of a heat bath and describes the 
interconnectivity between the cells. The higher h is set the lower is the in- 
terconnectivity of the network. The first term of Lp reflects the flip process 
Ui — — 1 — > Gi — +1 whereas the second term represents = +1 — ► <y% — — 1. 
The second part, Le-, expresses the exchange process Vi + Vj <-> Vj + Vi related 
to the existence of an active link between adjacent cells. The term in Eq.(||) 

A mi ■ ■ ■ A mf (10) 

<m\---m,f,i> 

represents the kinetic restriction which are mentioned above. The abbrevia- 
tion < mi ■ ■ ■ rrif,i > denotes the sets of all the / lattice cells neighbored to 
the cell i. The operator A m yields a non-zero result only if the cell m is an 
active one, so that the expression (|T0D differs from a zero value if it is applied 
to a cell of interest surrounded by at least / active cells. 
Using Eq.(^) the temporal evolutions of the two relevant observables (A^) 
and (Vfc) result in 

dt{A k ) = X B a X! {BkA mi ■ ■ ■ A mf ) — Xab X! (A k A mi ■ ■ ■ A m} ) 

<m\---nif ,k> <mi---rrif,k> 

d t (V k ) = 2D(Vk(A 2 k VkVk)) (11) 
where we exploit the discrete form of the Laplacian 

A k O k = ^J2(°r-O k ). (12) 

r(k) 
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The diffusion coefficient is modified, D = DqI 2 , where I is the length of the 
lattice cell. The current for the diffusive motion of the packages is given by 

j = -2DA 2 k VkV k - (13) 

It is intuitively obvious that the effective diffusion coefficient should be de- 
pendent on the squared concentration of the active cells, and therefore on 
the number of active links. 



3 Mean-field solution 

Because the formula for the vacancies concerns the same lattice index we 
may neglect it. Within the next step we apply a mean-field approximation 
and decouple all average values. Then, we obtain for the evolution equation 
of the packages 

d t (V(t)) = 2Ds7[(A(t))'s7(V(t))} (14) 
whereas the equation of motion for the mobile/active state yields 

d t (A (t)> = \ BA ((1 - A (t))(A (t)y - \ AB ((A (t))' +1 . (15) 
Its temporal solution is easily given by 

{A(t)) = A+[A(0)-A]exp^-^j (16) 



with the initial value A (0) and the steady state solution 

t Aba 1 



(17) 



Aab + Aba exp(2h) + 1 
As the inverse relaxation time of the flip process we find 

rf 1 = (Aba + \ab) CA f (18) 

with £ = z •...•(/ + 1). Notice that the steady state solution in mean-field 
approximation is equal to the solution for the paramagnetic lattice gas and 
independent of /. In contrast, the relaxation time depends on AK Because 
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the dynamics of the packages is globally conserved the steady solution is fixed 
for all time by the initial distribution, i.e. 



Testing the stability we analyze the Eqs.(|T4|) and (JfJ) in a linear manner. 
Assuming 5 A and 5V are small fluctuations around the steady-state values 
we obtain by means of (A (q,t)) = A + 5 A (q, t) and (V [q,t)) = V+ 5V (q, t) 
(q is the Fourier transformation of the local variable x) 

( 5A(q,t) \__( rf 1 \( 5A(q,t) \ 
0t [SV(q,t) )- ^ 2DA 2 q 2 ) { 8V (q,t) ) { > 

The divergence of the second relaxation time r^ 1 = 2DA 2 q 2 at the wave 
number q = reflects the global conservation of the information in our 
network. Due to the quadratic dependence on A the relaxation time T2 
rapidly increases if the active cells become passive. Obviously, the steady 
state reveals as stable against pertubations indicated by the negative sign. 
To gain more insight into the diffusion process of the information associated 



with packages we consider its evolution equation (II) applying the mean-field 
solution for the active state dl6|). This leads to a spatial independent, but 
temporal dependent effective diffusion coefficient with 



dt(V{t)) = 2D 



A+(A (0) - A) exp ' ' 



V 2 (V(t)). (21) 



A measure for the fluctuations of the packages concentration is given by 
F(t) = f\A(f)} 2 dt' (22) 







A 2 t + 2nA (A (0) - A) (l - e'n } + T ± [a (0) - A) 2 ( 1 



_2t_ 

e T i 

Studying the asymptotic limits of F (t) we recognize different temporal regimes. 
Whereas for small times the fluctuations are dominated by the initial value of 
the active cells F (t) ~ A 2 (0) t the fluctuations are approximated by F (t) ~ 
A 2 t for long times. Although the fluctuations go to infinity in the long-time 
limit (in this mean- field theory) there is a (quasi-) localization in an inter- 
mediate temporal range especially for small interconnectivity (equivalent to 
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low temperature). Only a small part of all cells is active (A <C (A (0) — A J) 
inhibiting the diffusion process. The exponential functions of the second and 
the third term in Eq. (^) are negligible so that the fluctuations F (t) remains 
almost constant. To see if there is a quasi-localization one must compare the 
relaxation time r in Eq. (|18|) and the time scale where the first term is equal 
the second and the third term in Eq. (|22|) . In this connection, we assume that 
the exponential functions may be neglected. The time £x associated with this 
point is (supposing A (0) > A) 



-2 




(23) 



Obviously, a localization it L ri) becomes only possible if the concentration 
A is small enough (sufficient low interconnectivity/temperature) for the time 
interval T\ < t < t^. A rough estimation yields the ratio 

— ~ exp (Ah) . (24) 

Further, we may use F (t) to transform the partial temporal derivative of the 
evolution equation, 

di ~ ITdFjf) ~ {A{t)) am' { ] 

and hence to obtain the ordinary diffusion equation 

d(V{t)) 



dF(t) 

applying the transformed Eq. (|H|) 



2-D V (V(t)) 



9 -^f L = ( A (t)) 2 %!§- = 2D(A(t)r V 2 (V(t)). (26) 

If we now start with a ^-distributed information (packages) density, i.e. we 
drop down the information at one point and observe how the information 
will be distributed in the time, 

V(x, 0) = V5 (x) (27) 
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we get the result of the Eq.flTJj) 



v & *) = : — ~vzj exp ■ ^ 

(4irDF(t))* \4xDF{t)J 

Hence, the behavior of F(t) directly influences the diffusion of the informa- 
tion. If F(t) remains almost constant in a time regime the diffusion process 
and therefore the distribution of the information stops, and data packages 
are localized (compare Fig.jl]). 




Figure 1: The diffusion of the packages concentration in the restricted net- 
work model (upper "localized" graph) in comparison to the ordinary diffu- 
sion process (lower graph) embedded in one dimension (D — 1, / = 1 and 
A(0) = |). 



4 Conclusions 

In our extended facilitated kinetic Ising model we find a (quasi-) localization 
due to the coupling of the diffusive dynamics for data packages with the 
existence of the active cells at sufficiently low temperatures. Because their 
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part of all cells is small in this case so that the diffusion almost stops. There- 
fore, the packages are fixed at or near the initial position, information cannot 
be widespread. Thus, we may separate two well-distinguished time regimes. 
First of all the fast flip (or active-passive) relaxation takes place. The higher 
the restriction number / the more the relaxation time slows down (compare 
Eq. (|T8|) ) . Then after a while of the relaxed system of active-passive states 
sticking the packages the influence of the (slow) diffusion is remarkable, equi- 
librates the packages concentration (information) of the total system. How- 
ever, for high enough temperature there are a sufficient part of active cells so 
that the coupled diffusion can fasten the equilibration and information can 
easily spread through the network. Like in |24| for the Fredrickson- Andersen 
model shown mean-field approximation may provide false results. The mean- 
field solution leads to a dynamics which completely breaks down below a 
critical (interconnectivity) temperature. That means below a critical (inter- 
connectivity) temperature information would spread through a fixed network 
like through a sponge (Bond percolation results would apply to this case, see 



e.g. Pq|.). Here, we exploit a more sophisticate approximation, taking local 
processes into account. But we expect that the found quasi- localization al- 
ters in a permanent localization at low temperature if / > z. Because thus 
it exist stable clusters of passive cells at any temperature in contrast to the 



case / < z where all cells can be activated EM. To prove this result one 



may study the theorem of van Enter describing a bootstrap percolation 
model |27j (or diffusion percolation |28 
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